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We discuss the Kretschmann, Chern-Pontryagin and Euler invariants among the second 
order scalar invariants of the Riemann tensor in any spacetime in the Newman-Penrose 
formalism and in the framework of gravitoelectromagnctism, using the Kerr-Newman 
geometry as an example. An analogy with electromagnetic invariants leads to the defi- 
nition of regions of gravitoelectric or gravitomagnetic dominance. 



1. Introduction 

This article shows how the most interesting second order scalar invariants of the 
Riemann tensor (i.e. the Kretschmann, Chern-Pontryagin and Euler invariants) 
among the set of 14 indepe ndent curvature scalar invariants can be expressed simply 
in any spacetime I ^ I ^ ISI not only in the Newman-Penrose formalism (NP), as is 
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well known in the literature, but also in the framework of gravitoelectromagnetism 
(GEM), the 1 + 3 splitting of spacetime. On the other hand, working instead with 
standard tensor algebra, the evaluation of such invariants involves long calculations 
for which a computer algebra system is often necessary.^ This widely held opinion 
is also reflected in a recent paper by Henry,^ which led to the present discussion 
and motivated further study. 

However, within the NP formalism, evaluating the second order scalar invariants 
of the Riemann tensor is not so difficult. Furthermore in special cases (Petrov type 
D spacetimes for example), this is almost trivial because the symmetries of the 
spacetime can be adapted to the NP frame (i.e. the two repeated principal null 
directions can be aligned with the two real null vectors of the NP tetrad) , with the 
subsequent vanishing of most of the Riemann tensor components. 

Similarly, adopting a GEM approach through a relative observer splitting of 
spacetime, the task of evaluating second order scalar invariants of the Riemann 
tensor can also be accomplished in a simple and elegant way which also allows 
their interpretation in complete analogy with the well known second order scalar 
invariants of the Maxwell 2-form. This is an an original contribution of the present 
article leading to a new interpretation of the regions in which such "quadratic 
curvatures" become negative. 

The following second order scalar invariants of the Riemann tensor will be stud- 
ied: the Kretschmann invariant ,(K i ) , and its "dual" counterparts which were omit- 
ted in the discussion of HenryM-lfhe Chern-Pontryagin (-K2) and the Euler invari- 
ant (# 3 )P These are discussed in detail for a Kerr-Newman spacetime, which is 
presently of interest in view of its possible astrophysical applications.^ 

Very recently such scalars have been considered in the context of black hole 
collisions and perturbations by Baker and Campanelli,!'-' who introduced a 'spe- 
ciality index' (involving second and third order invariants) which provides a mea- 
sure of how distorted a black hole perturbation is from the background space- 
time. During recent decades curvature invariants have also been considered in 
the context of quantum gravity and for effective theories of gravity applied to 
cosmology. In such cases, th ey play an important role: a) for the one-loop level 
rcnormalization of gravity^^ and b) for producing inflationary behaviour in early 
universe cosmology. Finally they have also been studied in the context of con- 
formal gravity, where in the post-Newtonian approximation limit, they generate 
corrections to the standard Newtonian potential which might have astrophysical 
implications}^!^ as a consequence of certain new gravitationa l fie ld equations in 
which the metric is coupled to matter fields via the Bach tensor.^l 

2. Kretschmann, Chern-Pontryagin and Euler Invariants 

We will consider the following second order scalar invariants of the Riemann tensor: 

Ki — R a p 7 sR af3jS , 
K 2 = [*R] aMS R aM \ 
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K 3 = [*R*U lg R af315 . (1) 

Starting from the relation between the Riemann, Weyl and Ricci tensors (or 
first Matte-decomposition of the Weyl tensorPD) 

RaPyS = CafiyS + ^(dayRpS ~ 9l3fRaS ~ OaSR^f + 9psR a j) 

-^(9^9/38 - g a 59p-i)R , (2) 
one can write K\ , K% and K 3 in the form: 

K x = C al)lS C a ^ s + 2R a pR a/3 - ii? 2 

= h + 2R a pR al3 - \r 2 , (3) 

= h , (4) 
K 3 = ~C aPl sC a ^ 5 + 2R a0 R al3 - ^R 2 

= -h + 2R a pR af3 -^R 2 . (5) 

where I\ and I2 are the only independent second order scalar invariants of the Weyl 
tensor: 

h = C aPl& C aP ~< & , (6) 
h = [*CU l5 C a ^ s , (7) 

because of the self-duality property of the Weyl tensor: [*C*] = —C. The relation 
K2 = I2 can easily be derived by using the tracefree property of the Weyl tensor 
and the Bianchi identities of the first kind for the Riemann tensor. It is worthwhile 
noting the following properties of K\, K2 and K 3 and I\ and I2: 

(1) The sum 

K x + K 3 = AR aP R afs -R 2 = An 2 [T a0 T afs -\t 2 ] 1 (8) 

as a consequence of the Einstein equations, cleanly separates off the non-Weyl 
part of the curvature determined by the matter content. Analogously one has: 

K x - K 3 = 2h + ^R 2 . (9) 

(2) The density associated with K 3l i.e. /C3 = ^/—gKs, is proportional to the 
topological Euler density: 

1 „ 1 



£ = "12^ K3 = J^V^iR^sR ^ - 4R a pR aP + R 2 ) ■ (10) 
Furthermore /C3 is the divergence of a vector density: 

£3 = -d a V a , (11) 
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where in coordinate components one has 



V a = 



-9V 



.1 



1 



(12) 



L 2 " ,u 3 

(3) if 2 & n d ^3 are two topological invariants obtained (in four dimensions) from 
the curvature two-forms. The Gauss-Bonnet theorem states that the integral 
of Kg over a compact manifold without boundary is proportional to its Euler 
characteristic. The integral of if 2 instead is related to the so called "instanton 
number" of the manifold!"! 



(4) The variational derivative of the density associated with I\, i.e. T\ 
is the Bach tensoi^l 



D 



a/3 _ 



1 51 X 



(13) 

V ~9 °9<x0 

which is symmetric and tracefree. The previous considerations concerning the 
topolo gica l invariants allow a simple and elegant way to calculate the Bach 
tensor (121 In fact from Eq. (gj it follows that 

i, = V^g(R af ^sR a0lS - m aP R af} + \r 2 ) , (u) 



which is equivalent to 

X x = 128tt 2 £ 



rl g(2R a pR 



a/3 



(15) 



The first term in l|15|) is the Euler density (|10fl . and consequently it does not 
contribute to Eq. G3J. 18 The variation of the second term gives directly the 
Bach tensor: 

2 11 

B a (3 = - R;a;0 ~ 2 R a 0;^'^ + — 9af3 R;ti ~ ^ R 9aj3 
4 



R Ra(3 + 9a(3 R)lv R^ ~ 4 R a ^S3v 



(16) 



2.1. Newman- Penrose Formalism 

The two second order scalar invariants of the Weyl tensor ii and I2 satisfy the NP 



relation 
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1 = 1,-112 = 16(3^| + ^4 " 4^3), 



(17) 



where for the NP notation we follow the conventions of ChandrasckharK^By taking 
the real and the imaginary parts (the latter reversed in sign) of Eq. I|17|l one finds: 



h = $e(I) 



(3</| + ip ip 4 - 4^3) + (3^1 + - H>iih) 



h = -9fm(J) = 8z (toft + VoV>4 - 4^ 3 ) - (3^1 + ^oV>4 - 4^ 3 ) • (1? 



Projecting R a pR a " onto the NP frame gives: 

RapR a ^ — 2[ i?n i?22 + -Rl2 2 — 2 i?i 3 i?24 ~ 2 i?i4 i?23 + i?33 R 



R, 



34 



= 8[ $00*22 + 18 A 2 + 2 $n 2 - 2 $oi$2i - 2 $10*12 + $02*20] ■ (19) 
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Finally R in the NP formalism becomes: 

R = 24A. (20) 

By inserting JTSJ, (|T§j) and lf2*UI) into and ©, we can express K\, K 2 and 

in terms of the NP quantities in a form valid in any spacetime. 

2.2. Gravitoelectromagnetic Formalism 

A 1 + 3 splitting of the spacetime is accomplished (locally) by introducing a family 
of test observers, i.e. a congruence of timelike lines with unit tangent vector u 
(u ■ u = —1). u defines both a local time direction and a local space (through its 
orthogonal subspace in the tangent space). Projecting tensor and tensor equations 
along u or orthogonally to u defines a 'measurement process' associated with the 
observer family it; for example the 'measurement' of a (j) tensor S results in a set 
of four spatial fields (i.e. for which any contraction by u vanishes): 

{u\S~< s , P(u) a lU 5 S~<s, P(u) 5 a u^ s ,P(u)%P(u) 5 ^ s }, (21) 

where P{u) a p = 5 a p + u a up is the projector orthogonal to u. 

Details about the splitting process systematically applied to spacetime tensor 
and spacetime differential operators as well as the origins of gravitoelectromag- 
netism can be found in Bini et al^^to which we refer for notation and conventions. 
Here we are interested in the splitting of the Riemann tensor. The only nonvanish- 
ing spatial fields associated with the 'measurement' of R by the observer family u 
are the following^ 

£{u) a /3 = Ra^vU^u", 

F(u)a/3 = \v(u) Pa a V(ur U 0R p a^ ; (22) 

£(u), Ti(u) and T{u) have been respectively called the electric-electric, the electric- 
magnetic and the magnetic-magnetic parts of the Riemann tensor. A straightfor- 
ward calculation shows: 

Kx = 4[Tr£(u) 2 - 2TrW(u) ■ H(uf + Tr^(u) 2 ] , 
K 2 = -8[TrW(u) ■ (£(u) - ^(u))] , 

K 3 = 8[Tr£(u) • F(u) + TrH(u) 2 ] . (23) 
3. Application: The Kerr-Newman Spacetime 

The Kerr-Newman solution is of Petrov type D and in an NP frame adapted to the 
two repeated principal null directions (a Kinnersley tetrad 12*1) one has only ip2 7^ 0, 
so that eqs. (|18|l become: 

h = 24(^ 2 2 h = 24i(^ 2 - ) . (24) 
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Moreover, as an electrovac solution of the Einstein-Maxwell system, the Kerr- 
Newman solution has a tracefree Ricci tensor R = as a consequence of the trace- 
free property of the energy-momentum tensor of the electromagnetic source field. 
Thus, in order to evaluate K\, K 2 and K 3 we need only calculate R a pR al3 or the 
components of the Ricci tensor, namely A and $ mrl . However, here A = and 
$ mn = 20 m n , the only surviving component of which is $n. Thus 

$n = -\{Ri2 +R34), A = - — (R12 - R34) = 0, (25) 
which can be inverted yclding 

R12 = R34 = -2(*u) = -4&?T. (26) 
By using these relations in eq. I|19|) . we have: 

R al3 R afi = 640 2 ^ 2 (27) 

so that 

K x = 24(tJj 2 2 + V|) + 128(0i^) 2 . (28) 

where P3I 

9= - (r -ia cos 6)-\ V2 = p 3 (M + Q 2 p) , 0i = ^Qp 2 . (29) 



The final result is: 



Ra0R «P = A ( 2^ T 2flH ( 30 ) 

(r 2 + cos^ 0) 4 



and 



Ki = 



(r 2 + a 2 cos 2 8) 
x [6A/ 2 ( r 6 - 15 r 4 a 2 cos 2 (9 + 15 rV cos 4 8 - a 6 cos 6 0) 
- 12MQ 2 r (r 4 - 10 r 2 a 2 cos 2 + 5a 4 cos 4 0) 

+ Q 4 (7 r 4 - 34 r 2 a 2 cos 2 6 + 7 a 4 cos 4 0) ] , (31 ) 

which coincides with the expression obtained by Henrjd using a computer algebra 
system. In the same way we can calculate K2 = 24i(^f — ipV)'- 

96a cos 



Ki 



(r 2 + a 2 cos 2 0) 
x [M 2 r (3r 4 - 10 r 2 a 2 cos 2 8 + 3a 4 cos 4 1 



MQ 2 (5 r 4 - lOaV cos 2 6 + a 4 cos 4 0) 

2rQ 4 (r 2 - a 2 cos 2 0)] , (32) 
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and K 3 = -24(^1 + + 128(</>i0i) 2 : 



if a = 



(7- 2 + a 2 cos 2 0) 6 



[6M 2 ( r 6 - 15 r 4 a 2 cos 2 + 15 rV cos 4 (9 - a 6 cos 6 6) 
12MQ 2 r (r 4 - 10 rV cos 2 (9 + 5a 4 cos 4 6) 

g 4 (5 r 4 - 38 rV cos 2 (9 + 5 a 4 cos 4 0) }. (33) 



The expression of 1\ (I 2 = K 2 ) is instead 
48 

(r 2 + a 2 cos 



48 

J i = 7Z5 , „2„„ 2 m6 ( r4 + » 4 cos4 ~ 6fl2 cos2 ^)(Q 2 -Q+)(Q 2 -Q-) (34) 



iC 3 = -48 + 96 ™# - 8 - 126 ° + Q ( \ 

Instead, expanding for small a/M and Q/M one finds 



(35) 



with 

^ 2 M(r + acos6>)(r 2 + a 2 cos 2 ± 4arcos6>) 
± r 2 _ a 2 cog 2 $ _|_ 2 ar cos $ 

At large distances K\, K2 and K3 have the expansion!^! 
gl = 48 - 96 g£ + 56 ^ ~ 18 y Q2 C ° S2 * + Q (1 

K 2 = 288 M \ CQS * - 480 MQ ™ 6Q2 + l92 aC ° S9 & = " ^ 3 + Q ^ 1 



in 



Af 2 M 2 a 2 cos 2 AfQ 2 
Ki = 48 — - 1008 5 96 — + ... ~ -Jf 8 

M 2 a cos Af 2 a 3 cos 3 Ma cos 0Q 2 
#2 = 288 = 2688 5 480 =— ^- + ... 

It is also possible to find the surfaces on which K\, K 2 and K3 vanish. Solving 
K 2 = for cos0 in the Kerr-Newman spacetime gives the following curves: 

cos 6 = 0, 

y/Mr(3Mr - 2Q 2 ) 



cos = ±- 



Ma 



C ° S9 = ± a(3M,-Q 2 ) ' (36) 



a According to the notation of the present paper, the quantity K introduced by de Felice as 
the "first Khretshmann invariant" actually is I\. 
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which, in the Kerr case, reduce to: 

cos = 0, 

cos0 = ± r, 

a 

cos9 = ±—^r, (37) 

To solve K\ = it is convenient to re-express it in terms of the new variable 
Z = a 2 cos 2 0, yelding the cubic equation: 

-6 M 2 Z 3 + (90 r 2 M 2 - 60 MrQ 2 + 7 Q 4 )Z 2 



(r 2 + Zf 1 

+ (-90 r 4 M 2 + 120 Mr 3 Q 2 - Mr 2 Q 4 )Z + 6M 2 r 6 

-12 Mr 5 Q 2 + 7 Q 4 r 4 ] = , (38) 

which can be solved exactly by using the Cardano formulas. Instead of giving these 
lengthy formulas here, we give the simpler solutions for the Kerr case: 

r = ±acos0, 

r = (2 ± V3) a cosy, 

r = (-2±V3)acos6 . (39) 

The equation K% = may be handled similarly; for the simpler Kerr case it reduces 
to K3 = —K\ = and has the same solutions. 

Figs. JTJ-@I and ©~© show these surfaces for the Kerr and Kerr-Newman 
cases respectively for two typical parameter choices. As usual in the literature 
the plots have been drawn in polar-like coordinates (x,z), related to the Boyer- 
Lindquist coordinates (r, 9) by the relations x = r sin 9 (horizontal axis), z = r cos 9 
(vertical axis). Note that in the (x, z) coordinates, curves represented by the (polar) 
relation r = Acos#, (with A constant) as in Eqs. (|39|l correspond to circles with 
center at (0, A/2) and radius A/2. 

We conclude this section by listing the gravitoelectromagnetic components of 
the Riemann tensor in the Kerr-Newman spacetime, which once inserted into eqs. 
Ij23(l allow an alternative (to NP) and straightforward derivation of the second 
order scalar invariants of the Riemann tensor. The Carter observer family with 
four-velocity 

(r 2 + a 2 ) 



u. 



(car) = 7= (dt + "J— 2^) , (40) 



'AS 

as expressed in Boyer-Lindquist coordinates, where A = r 2 + a 2 — 2Mr + Q 2 and 
S = r 2 + a 2 cos 2 9, induces the 1 + 3 decomposition of the orthonormal frame natu- 
rally associated with the Kinnersley null frame. One finds the following coordinate 
components with respect to the Boyer-Lindquist coordinates 

n 2 sin 2 9 

^ («(car))« = ^r^[ a2 (« 2 C0S2 0) + Mr ( r2 ~ 3fl2 C0S2 )1 
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— ~ 2^( U (car))t?i 



a 2 sin 2 i 
£ 2 



■£(w(car))e0 
7^( M (car))^! 



(r 2 + a 2 ) 2 

2Mr(r 2 - 3a 2 cos 2 0) + Q 2 {a 2 cos 6 -3r 2 ) 

£ [U( caI ))rr = ^2 ' ( 41 ) 



and 

2Q 2 

•^("(car))^ = -£("(car))a/3 + j^^atf} , 

v _ a cos fl[.M(a 2 cos 2 6 — 3r 2 ) + 2rQ 2 ] Q 2 

H( U(car) ) Q/3 - - Q 2(a 2 C0S 2 ,2) + Mr(r 2 _ 3fl2 C0S 2 fl) ^("WW + YA. 13 



4. Gravitoelectric or Gravitomagnetic Dominance in Vacuum 
Spacetimes 

For vacuum spacetimes the Riemann and the Weyl tensors coincide and the relation 
K\ = —K3 holds. As a consequence in a 1 + 3 spacetime splitting point of view, 
the magnetic-magnetic part of the Riemann tensor reduces to its electric-electric 
part (apart from a sign, i.e. J~{u) = —£{u)), while the electric-magnetic part Ti(u) 
becomes symmetric. This is a simplified situation in which one has only two in- 
dependent spatial tensor fields representing the Riemann tensor, exactly as in the 
electromagnetic case where the Maxwell 2-form F is represented by the electric and 
magnetic vector fields; thus it is possible here to push forward the analogy between 
electromagnetism and gravitoelectromagnetism by considering the correspondence 
between the invariants of F and those of R. 

The detailed analysis of the Maxwell invariants in a curved background dates 
back t o th e seventies and to the pioneering work of Ruffini, Hanni, Damour and 
WilsonPBl They introduced the concepts of a 'region of electric dominance' (where 
for any observer u, \\E(u)\\ > ||S(u)||, i.e. F^F^ < 0), that of 'region of magnetic 
dominance' (where for any observer u, \\E(u)\\ < \\B(u)\\, i.e. F^ U F^ U > 0) and 
that (stronger than the first two) of a 'plasma horizon,' i.e. the boundary of a region 
in which the magnetic field can support an infinitely thin plasma against Coulomb 
attraction. In this case the meaning of the various regions in terms of trapping of 
particles by the magnetic field is clear. 

When working out the analogy with the invariants of the Riemann tensor the 
most natural candidate to play the role of F^F^ is K\. One can rephrase the 
Ruffini et al definitions in terms of the corresponding gravitoelectromagnetic quan- 
tities by introducing regions of 'gravitoelectric dominance': Tr£(u) 2 > Tr7i(u) 2 , 
i.e. K\ > 0, and regions of 'gravitomagnetic dominance': Tr£(u) 2 < Tr7i(u) 2 , i.e. 
K\ < 0. This gives a more reasonable interpretation of the regions where such a 
'quadratic' curvature becomes negative rather that claiming it to be 'a new type of 
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curvature' as does Henry™ 

In Figs. (l)-(4) the regions of gravitoelectric and gravitomagnetic dominance 
can be easily recognized. To this end one must consider the dashed circles which 
correspond to the vanishing of K%: as a general feature, far from the hole the 
spacetime is gravitoelectrically dominated (K\ > 0) while a gravitomagnetically 
dominated region (K\ < 0) exists close to outer horizon (in other words, K\ is 
positive at the spatial infinity and changes sign at each crossing of the dashed 
circles). 

Moreover, our discussion is in agreement with a pioneering work of de FelicJ^ 
which introduces the concept of 'repulsive domains' in a curved spacetime (related 
to the definition of a 'negative effective mass' for the black hole) as possible markers 
of the existence of a timelike singularity in terms of the vanishing of K\ on surfaces 
containing the singularity itself. 

Finally the analog of the concept of plasma horizon must be substantially revised 
since Riemann tensor tidal forces act only on particles with intrinsic structure, like 
spinning test particles. In principle it is possible to attempt such a generalization 
but one must be careful in specifying what kind of particle/fluid trapping is under 
consideration. 

5. Conclusions 

The Kretschmann, Chern-Pontryagin and Euler invariants of the Riemann tensor 
have been written both in the NP and GEM formalisms for any spacetime in a form 
for which their evaluation can be done easily without the use of a computer algebra 
system, despite the widespread opinion in the literature to the contrary. These 
scalars have been examined in the Kerr and Kerr-Newman cases. The regions of 
spacetime where such objects vanish have been studied analytically, leading to the 
definitions of regions of gravitoelectric and gravitomagnetic dominance, at least for 
a vacuum spacetime. This approach is in complete analogy with the electromagnetic 
case and suggests an alternative interpretation of the scalar invariants. 
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Fig. 1. Kerr Black Hole with parameters M = 1, a = 0.5. Dashed lines in Fig. la correspond 
to both K\ = and K3 = (K\ = —K3 in vacuum); thin solid lines in Fig. lb correspond to 
K2 = 0; thick solid lines in both Fig. la and Fig. lb are the inner/outer horizons. 



M=l, a=0.5, Q=0 




Fig. 2. Kerr Black Hole with parameters M = 1, a = 0.5. Superposition of Figs, la and 
Dashed lines correspond to both K\ = and K% = (Ki = —K3 in vacuum); thin solid li 
correspond to K2 = 0; thick solid lines are the inner/outer horizons. 
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Fig. 3. Kerr Black Hole with parameters M = 1, a = 0.95. Dashed lines in Fig. 3a correspond 
to both K\ = and K3 = (K\ = —K3 in vacuum); thin solid lines in Fig. 3b correspond to 
K2 = 0; thick solid lines in both Figs. 3a and 3b are the inner /outer horizons. 
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Fig. 4. Kerr Black Hole with parameters M = 1, a = 0.95. Superposition of Figs. 3a and 3b. 
Dashed lines correspond to both K\ = and K3 = (Ki = —K3 in vacuum); thin solid lines 
correspond to K2 = 0; thick solid lines are the inner/outer horizons. 
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Fig. 5. Kcrr-Ncwman Black Hole with parameters M = 1, a = 0.6, Q = 0.6. Dashed lines in 
Fig. 5a correspond to K\ = 0, thin solid lines in Fig. 5b to K2 = 0, grey solid lines in Fig. 5c to 
K3 = and thick solid lines are the inner/outer horizons. 




Fig. 6. Kerr-Newman Black Hole with parameters M = 1, a = 0.6, Q = 0.6. Superposition of 
Figs. 5a, 5b and 5c. Dashed lines correspond to K\ = 0, thin solid lines to K2 = 0, grey solid lines 
to Kz = and thick solid lines are the inner/outer horizons. 
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Fig. 7. Kcrr-Ncwman Black Hole with parameters M = 1, a = 0.7, Q = 0.7. Dashed lines in 
Fig. 7a correspond to K\ = 0, thin solid lines in Fig. 7b to K2 = 0, grey solid lines in Fig. 7c to 
K3 = and thick solid lines are the inner/outer horizons. 
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M=l, a=0.7, Q=0.7 
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Fig. 8. Kerr-Newman Black Hole with parameters M = 1, a = 0.7, Q = 0.7. Superposition of 
Figs. 7a, 7b and 7c. Dashed lines correspond to K\ = 0, thin solid lines to K2 = 0, grey solid lines 
to Kz = and thick solid lines are the inner/outer horizons. 



